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Performing molecular dynamics simulations for all-atom models, we characterize the conforma- 
tional and structural relaxations of poly(ethylene oxide) and poly(propylene oxide) melts. The 
temperature dependence of these relaxation processes deviates from an Arrhenius law for both poly- 
mers. We demonstrate that mode-coupling theory captures some aspects of the glassy slowdown, 
but it does not enable a complete explanation of the dynamical behavior. When the temperature 
is decreased, spatially heterogeneous and cooperative translational dynamics are found to become 
more important for the structural relaxation. Moreover, the transitions between the conformational 
states cease to obey Poisson statistics. In particular, we show that, at sufficiently low temperatures, 
correlated forward-backward motion is an important aspect of the conformational relaxation, lead- 
ing to strongly nonexponential distributions for the waiting times of the dihedrals in the various 
conformational states. 



I. INTRODUCTION 

Understanding the glass transition of polymer melts is 
of enormous interest from the viewpoints of fundamental 
and applied science. When polymer melts or other glass- 
forming liquids are cooled, the structural relaxation slows 
down tremendously and, eventually, an amorphous solid, 
a glass, is formed at the glass transition temperature T g . 
The structural relaxation, or, equivalently, the a relax- 
ation of most glass-forming liquids exhibits two striking 
features: its time dependence differs from a single expo- 
nential function and its temperature dependence deviates 
from an Arrhenius lawi 

Molecular dynamics (MD) simulations have proven a 
powerful tool to investigate the initial stages of the glassy 
slowdown of the structural relaxation in moderately vis- 
cous liquids ,2iMi5i& Most of these studies focused on sim- 
ple model-glass formers, e.g., binary Lennard- Jones mix- 
tures. In MD simulations of glass- forming polymer melts, 
coarse-grained models were used, while investigations on 
chemically realistic models are more rare.^ 6 - For the latter 
models, consideration of all the interactions associated 
with the connectivity, e.g., of the energy barriers against 
the conformational dynamics, interferes with the neces- 
sity to follow the slowdown of the molecular dynamics 
over broad temperature and time ranges due to the lim- 
ited computer power. However, the rapid development of 
computer technology steadily improves the possibilities 
and, very recently, chemically realistic polymer models 
have started to become an important tool to study the 
mechanisms for the primary and Johari-Goldstein sec- 
ondary relaxation processes in polymer meltsi^ 

The mode-coupling theory (MCT), which focuses on 
density fluctuations, has been put forward to explain the 
observation that the temperature dependent structural 
relaxation times deviate from an Arrhenius law for most 
glass-forming liquids^ In its idealized version, MCT pre- 
dicts a power-law divergence of the a-relaxation time at a 



critical temperature T c . MD simulations were employed 
to test the predictions of MCT for moderately viscous liq- 
uids. It was found that this theory captures many aspects 
of the glassy slowdown for several simple model-glass for- 
mers, including a bead-spring polymer modele d ! 11 while 
the applicability to all-atom polymer models was contro- 
versially discussed£ & 12 ' 13 In any event, T c is substan- 
tially higher than T g , indicating that MCT fails to de- 
scribe the molecular dynamics in highly viscous liquids. 1 
Therefore, it is important to consider further aspects of 
glass-forming liquids and approaches focusing on the co- 
operativity and heterogeneity of the structural relaxation 
have received considerable attention . 14 ' 15 

Various experimental observations demonstrated that 
the nonexponential a relaxation of glass-forming liq- 
uids is related to the effect that the molecular dynam- 
ics are heterogeneous, i.e., it is possible to select par- 
ticles that rotate or translate much farther or shorter 
distances than an average particle. 16 ' 17 ! 18 However, most 
experimental techniques provide only limited information 
about the spatial distribution of particles showing differ- 
ent mobilities. Recent nuclear magnetic resonance ap- 
proaches demonstrated for various glass-forming liquids, 
including a polymer melt, that the dynamics are spa- 
tially heterogeneous . 19 ' 20 ' 21 Specifically, particles within 
a physical region of the liquid show an enhanced or dimin- 
ished mobility as compared to particles in a region a few 
nanometers away. Nevertheless, a detailed experimental 
characterization of the time and temperature dependence 
of spatially heterogenous dynamics is still lacking. 

MD simulations provide straightforward access to spa- 
tial correlations of the particle mobility^! 6 - Work on 
simple models of atomic and polymeric glass-forming 
liquids reported that highly mobile and highly immo- 
bile particles aggregate into clusters, which are tran- 
sient in natuTe^^^^L^l The clusters of highly 
mobile particles are largest in the very early stages of 
the a relaxation. Upon cooling, the cluster size in- 
creases and a divergence near the critical temperature T c 
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was proposecL ' 25 i 26 Moreover, it was demonstrated that 
string-like motion is an important channel for the struc- 
tural relaxation of the highly mobile particles . 29 ' 30 ! 31 
This means that mobile particles tend to follow each 
other along one-dimensional paths. The existence of spa- 
tially heterogeneous dynamics is not restricted to simple 
models, but this effect was also observed for more com- 
plex models of glass forming liquids, namely, for models 
of propylene carbonate, water, and silica. 32 ' 33 ' 34 ' 35 ' 36 In- 
terestingly, the cooperative string-like motion was found 
to be of little relevance for the silicon atoms in the silica 
model, where the structural relaxation follows an Arrhc- 
nius lawM& 

It was argued that the structural relaxation differs be- 
tween bead-spring and all-atom polymer models.— Specif- 
ically, conformational dynamics are of central importance 
for the latter, but not for the former models. For all- 
atom polymer models, the slowdown of the a relaxation 
was argued to depend not only on intermolecular pack- 
ing effects, as assumed in MCT, but also on intramolec- 
ular torsional barriers i 5 ' 8 ' 12 ' 13 Moreover, previous work 
showed that the conformational relaxation is by no means 
a simple relaxation process, but it involves, e.g., pro- 
nounced dynamical heterogeneities ^ 37 ' 38 In view of all 
these results, it becomes apparent that improving our 
understanding of the polymer glass transition requires a 
detailed characterization of the interplay of the intra- and 
intermolecular aspects of the molecular dynamics. 

Here, we perform MD simulations for all-atom models 
of poly(ethylene oxide) (PEO) and poly(propylene ox- 
ide) (PPO). Due to the capability to dissolve salts, PEO 
and PPO are popular materials for the preparation of 
polymer electrolytes^ Therefore, MD approaches stud- 
ied how the presence of ions affects the structure and 
the dynamics of these polymers both in the bulk and in 
confinement.^^^^^^^^iii&i 9 . Our investigation 
focuses on the temperature dependent dynamics of the 
neat model polymer melts. Specifically, we study the 
applicability of MCT and the relevance of spatially het- 
erogeneous and cooperative dynamics for the structural 
relaxation. Moreover, we perform a detailed characteri- 
zation of the conformational relaxation, elucidating the 
importance of correlated forward-backward jumps. 



£ V to (^ kl ) + V nb ({r}). (I) 

dihedrals 

Here, {r} is the set of all atomic coordinates. The bonded 
interactions are comprised of energies due to stretch- 
ing of bonds, bending of valence angles, and torsion of 
dihedral angles. The nonbonded interactions V nb are 
composed of Coulombic and van der Waals interactions, 
the latter being modeled using a Buckingham potential. 
The explicit form of the various interaction terms and 
the corresponding potential parameters are given in the 
literature. 50 ' 51 ' 52 In the case of PEO, we apply the force 
field termed FF-3 in Ref. It was demonstrated that 
these models enable a good reproduction of thermody- 
namic, structural, and dynamical aspects of PEO and 
PPO meltsM&dl 

The MD simulations were performed using the GRO- 
MACS software package.— We applied periodic bound- 
ary conditions and a time step of 1 fs. The nonbonded 
interactions were calculated utilizing a cutoff distance of 
12 A. To treat the Coulombic interactions, the particle- 
mesh Ewald technique was employed.— The LINCS al- 
gorithm was used to constrain all bonds ^ Prior to 
data acquisition, the systems were equilibrated in sim- 
ulations at constant N, P, and T, using the Rahman- 
Parrinello barostat 56 and the Nose- Hoover thermostat. 57 
These equilibration runs, which spanned 20—30 ns at the 
lower temperatures, allowed us to adjust the densities 
p(T). The density increases from p(450K) = 1.01 g/cm 3 
to p(280K) = 1.12 g/cm 3 for PEO and from p(450K) = 
0.89 g/cm 3 to p(250K) = 1.05 g/cm 3 for PPO. The sub- 
sequent production runs were performed in the canonical 
ensemble, i.e., at constant TV, V, and T, employing the 
Nose- Hoover thermostat. Although experimental work 
demonstrated that PEO is partially crystalline at room 
temperature and ambient pressure^ 9 , the simulation re- 
sults give no evidence for an onset of crystallization in 
the studied time and temperature ranges. 

III. RESULTS 

A. General characterization of the structural 
relaxation 



II. METHODS 

The studied PEO and PPO models are comprised of 
polymer chains, H-[CH 2 -0-CH 2 ]i 2 -H and CH 3 -0-[CH 2 - 
CH(CH3)-0]n-CH3, respectively. Each model is com- 
posed of 32 chains. The PPO chains are atactic, i.e., the 
side groups are randomly connected. The interatomic in- 
teractions of these polymer models are described by two 
well established quantum-chemistry based, all-atom force 
fields ; 50 ' 51 which can be written in the form 

V({r}) = £ V b °{r l3 )+ £ V»e(e tjk ) + 

bonds angles 



It is well established for the a relaxation of polymer 
melts that the time dependence differs from a simple ex- 
ponential behavior and the temperature dependence does 
not obey an Arrhenius law.— Therefore, we first demon- 
strate that the studied PEO and PPO models show these 
key features. While the incoherent intermediate scatter- 
ing function 

F s (q,t) = (cos{q - [rfo+t)-r(t )]» (2) 

provides us with information about translational motion, 
the orientational autocorrelation function 

F 2 (t) = ^(3[e(t + t)-e(t )] 2 -l) (3) 
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FIG. 1: Temperature dependent incoherent intermediate scat- 
tering functions F s (q, t) for the oxygen atoms of PPO and 
q = 0.91 A . The temperatures are, from left to right, 450 K, 
400 K, 350 K, 325 K, 300 K, 280 K, 265 K, and 250 K. The 
dashed lines are von Schweidler fits, see Eq. J9|). The inset 
shows the time temperature superposition of the same data. 
The time constants t from KWW interpolations of the a re- 
laxation were used for the scaling of the time axis. The dashed 
line is a KWW fit of the a relaxation ((3 = 0.72). 



intermediate scattering functions F s (q,t) for the oxygen 
atoms of PPO show a pronounced temperature depen- 
dence, in particular at low temperatures. To quantify the 
slowdown of the structural relaxation, we extract temper- 
ature dependent translational and rotational correlation 
times according to F s (q,TT) = 1/e and ^(tr) = 1/e. 
Figure shows the results for PEO and PPO. For both 
models, tt and exhibit a comparable temperature de- 
pendence that cannot be described by an Arrhenius law, 
as expected for polymer melts. Rather, a Vogel-Fulcher- 
Tammann (VFT) law^ 

T (T)=r 00 exp(-^j , (4) 

enables good interpolations of the data. VFT fits to tt 
(tr) yield T = 179 K (T = 191 K) for PEO and T = 
179 K (T = 175 K) for PPO. The results for PPO are in 
reasonable agreement with To = 157— 170 K obtained in 
light scattering and dielectric spectroscopy studies . 60 ! 61 

Inspecting the shape of the scattering functions in Fig. 
[U we see nonexponential decays. As expected for glass- 
forming liquids^ the time dependence in the a-relaxation 
regime is well described by a Kohlrausch- Williams- Watts 
(KWW) function, or, equivalently, stretched exponential 
function 



yields insights into rotational motion. Specifically, the 
scattering function F s depends on the translational dis- 
placements [r(io+i)— r(io)] of the atoms during the time 
interval t, where the absolute value of the scattering vec- 
tor, <jf=|q|, determines the length scale on which dynam- 
ics is probed. Throughout this contribution, we focus on 
the translational motion of the oxygen atoms, but we en- 
sured that qualitatively similar findings are obtained for 
the carbon atoms. The orientational correlation function 
F 2 depends on the angular displacements \e(t +t) -e(t )| 
during the time interval t. We study the reorientation of 
the C-H bonds and e(t) is the unit vector describing the 
direction of a C-H bond at time t. 2 H NMR stimulated- 
echo experiments were used to measure i 7 ^) for C-D 
bonds of deuterated PPO.— Mimicking the experimen- 
tal situation, we restrict the analysis to C-H bonds in 
the methylene groups to avoid effects from fast three- 
fold methyl group jumps. Finally, in Eqs. ([2]) and ([3]), 
the brackets (...) denote the average over various time 
origins to and over all atoms or bonds belonging to the 
considered atomic or bond species. 

First, we use these correlation functions to ascertain 
the temperature dependence of the translational and ro- 
tational motion associated with the a relaxation. To ad- 
dress the translational aspect, we calculate F s (q,t) us- 
ing q = 1.31 A" 1 and q = 0.91 A" 1 for PEO and PPO, 
respectively. These values of the momentum transfer 
correspond to the respective position of the first maxi- 
mum of the intermolecular oxygen-oxygen pair distribu- 
tion functions. 44 ' 45 In Fig. [TJ we see that the incoherent 



FKWwit) = A exp 



(0</3<l). (5) 



Here, r and (3 quantify the time scale and the stretching, 
respectively. Fitting F s (q, t) to a KWW function, we find 
stretching parameters /3~0.63 for PEO and /3~0.71 for 
PPO, independent of temperature in the studied tem- 
perature ranges. For PPO, the time-temperature super- 
position is further demonstrated in the inset of Fig. Q] 
Clearly, the curves coincide in the a-relaxation regime 
when the time axis is scaled with the correlation time 
t. KWW interpolations of T^i) yield stretching param- 
eters (3 fa 0.36 for PEO and (3 pa 0.48 for PPO, consis- 
tent with stretching parameters obtained from 2 H NMR 
stimulated-echo experiments near T g £^ Thus, with re- 
spect to both temperature and time dependence, the 
structural relaxation of the studied polymer models re- 
sembles that of PEO and PPO melts, confirming the 
quality of the used force fields 4 s - 

We note that, for PPO at T = 250 K, the transla- 
tional and rotational correlation functions give evidence 
for some deviations from time-temperature superposition 
in the early stages of the decay, where the curves de- 
crease significantly faster than expected from the KWW 
behavior, see Fig.[TJ Since dielectric spectroscopy studies 
demonstrated that PPO exhibits a Johari-Goldstein sec- 
ondary relaxation, which starts to separate from the pri- 
mary relaxation at T ps 250 K ) 60 i 61 i 62 we attribute these 
deviations to the onset of this secondary process. 
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FIG. 2: Time constants characterizing the dynamics of the 
(a) PPO and (b) PEO models. The non-Gaussian parameter 
Q2, the mean cluster size S w and the mean string length L w 
are a maximum at t a 2, ts and t^, respectively. Moreover, 
F a (q, tt) = l/e and ^(tr) = 1/e. The solid and dashed lines 
are VFT interpolations of r T (T) (PPO: T = 179 K, PEO: 
T = 179 K) and r R {T) (PPO: T = 175 K, PEO: T = 191 K). 



B. MCT analysis 

The time temperature superposition observed for the 
a relaxation of the PEO and PPO models is in agree- 
ment with MCT. 9 Therefore, we now check whether fur- 
ther predictions of this theory are fulfilled. According to 
MCT, the particles of glass-forming liquids are trapped 
in cages formed by their neighbors for some time until an 
escape from these cages is possible during the structural 
relaxation. As a consequence of this trapping, a plateau 
regime, or, equivalently, /3-relaxation regime, preceding 
the a-relaxation regime develops when the temperature 
is decreased towards the critical temperature T c . Quan- 
titatively, MCT predicts a power-law divergence of the 
ct-relaxation time at the critical temperature T c : 
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FIG. 3: Test of MCT predictions for PPO on the basis of 
the incoherent intermediate scattering functions F s (q, t) of 
the oxygen atoms, (a) Correlation times T q resulting from 
KWW fits of the a-relaxation regime for the indicated values 
of q. The solid lines are MCT power laws r q oc (T— T c ) ~ 7 where 
T c — 242 K and 7 = 2.8. (b) Temperature dependence of h^/ lh , 
using 76= 1.44 for g = 0.50 A -1 , 76= 1.34 for q = 0.91 A -1 , 
and 76 = 1.16 for q = 2.00 A -1 . The solid lines are linear 
interpolations h 1 q hb oc (T-242 K). 



Another key prediction of MCT is the factorization the- 
orem for the /3-relaxation regime. It states that, for such 
times, all correlation functions, in particular, the incoher- 
ent intermediate scattering functions for different values 
of the momentum transfer, can be written as 



F s (q,t) = f c q +a q G(t). 



(7) 



Here, the plateau value f° and the amplitude a q depend 
on the value of q, while the (3 correlator G(t) is indepen- 
dent of the observable. If the factorization theorem is 
obeyed, the curves 



R(q,t) = 



F s (q,t)-F s (q,t>) 
F s (q,t") - F s (q,t') 



(8) 



r(T) oc (T - T c )~\ 



(6) 



for different values of the momentum transfer will col- 
lapse onto a master curve, provided the times t' and t" 
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are chosen inside the /3-relaxation regime^, Hence, cal- 
culation of R(q, t) allows one to check the factorization 
theorem without invoking a fitting procedure. G(t) can 
be expanded for times close to the central /^-relaxation 
time t a . The expansion for t > t a leads to the von Schwei- 
dlcr law 



Fs(q,t) = f c q 



(9) 



showing that a power law characterized by the universal 
von Schweidler exponent b describes the initial stages of 
the decay from the plateau. For temperatures T > T c , 
the amplitude h q decreases upon cooling according to 



h q (T) oc (T 



76 



(10) 



Finally, within MCT, b and 7 are related via the exponent 
parameter A: 



A = 



r 2 (i- Q ) _ r 2 (i + &) 

r(l-2a) ~ T(l + 26) 
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To check these MCT predictions for the PPO model, 
we analyze the scattering functions F s (q,t) of the oxy- 
gen atoms for various values of the momentum transfer q. 
First, we fit a KWW function to the scattering functions 
in the a-relaxation regime. Figure [3ja) shows the tem- 
perature dependent correlation times r q resulting from 
these fits for three values of q. A MCT power law, see 
Eq. ©, with T C = 242±2K and 7 = 2.8±0.1 nicely de- 
scribes all data at T/T c — 1 > 0.1. These findings are in 
reasonable agrement with T C = 236±2K and 7 = 3.7±0.8 
from experimental studies i 60 ' 64 However, there are devia- 
tions from the MCT power law at the lowest temperature 
T/T c — 1«0.03. At the present, it is not clear whether 
these deviations are related to the possible onset of the 
Johari-Goldstein secondary relaxation near T c . Qualita- 
tively similar deviations from the MCT predictions were 
reported for a chemically realistic model of polybutadiene 
in the vicinity of the critical temperature^ 

Next, we fit the von Schweidler law to the scattering 
functions F s (q,t). Due to a possible interference of the 
Johari-Goldstein secondary relaxation, we exclude the 
data for the lowest temperature T = 250K from further 
MCT analysis. The von Schweidler law enables a good in- 
terpolation of the decays in the late-/3/ early-a relaxation 

regime, see Fig. H] For g = 0.91 A , the von Schweidler 
fits yield / g c = 0.93±0.01 and 6 = 0.46±0.01 independent of 
temperature. According to Eq. pip , this value of the von 
Schweidler exponent b translates into 7 = 2.9, in agree- 
ment with 7 = 2.8±0.1 determined from the temperature 
dependent correlation times r q . To check the validity of 
Eq. (p~0|) . we use the exponents & = 0.46 and 7 = 2.9 and 
plot the temperature dependence of hl^ jb in Fig. E^b). 
In harmony with the MCT prediction, there is a linear 
relationship and, by extrapolation, h,y^ b vanishes at the 
critical temperature T c = 242 K, extracted from the tem- 
perature dependence of the a relaxation. 
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FIG. 4: Test of the MCT factorization theorem for PPO at 
T = 265K. The functions R(q,t) were obtained from the data 
F B (q, t) for the oxygen atoms according to Eq. JSp. The values 
of the momentum transfer q are indicated and times t" = 2 ps 
and t' = 100 ps were used. 



However, the MCT predictions are not fulfilled when 
including the scattering functions for other values of q 
into the analysis. Von Schweidler fits for q = 0.5 A and 
g = 2.0 A" 1 yield 6 = 0.59±0.02 and 6 = 0.23±0.02, respec- 
tively, corresponding to 7 = 2.4 and 7 = 5.0. Thus, there 
is no universal von Schweidler exponent and, for high 
and low values of q, the values of 7 calculated from the 
von Schweidler exponents deviate from 7 = 2.8 resulting 
from the temperature dependence of the a relaxation. 
To demonstrate the violation of the factorization theo- 
rem independent of any fitting routine, we show R(q, t) 
for different values q in Fig. 2J The data was calculated 
from F s (q,t) for the oxygen atoms of PPO at T = 265K. 
Clearly, the data do not collapse onto a master curve in 
the /J-relaxation regime, indicating the violation of the 
factorization theorem. Qualitatively similar results were 
observed at all studied temperatures. 



C. Spatially heterogeneous dynamics 

Other models of the glass transition focus on the het- 
erogeneity and the cooperativity of the dynamicsj 14 ! 15 
The importance of these effects was demonstrated 
in MD simulation studies on various glass-forming 

liqttjds 22.23.24.25.26.27.28.29.30.31.32.33.34.35.36 JJ ere we m _ 

vestigate the relevance of heterogeneity and coopera- 
tivity for the first time for chemically realistic poly- 
mer models. For various model-glass formers, it was 
found that the distribution of scalar particle displace- 
ments \r(to+t)—r(to)\ deviates from a Gaussian at inter- 
mediate times t between ballistic and diffusive motion. 
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FIG. 5: Non-Gaussian parameter a,%(t) for the oxygen atoms 
ofPPO at various temperatures (450 K, 400 K, 350 K, 325 K, 
300 K, 280 K, 265 K, 250 K). 



These deviations, which are a first indication for the ex- 
istence of heterogeneous dynamics, can be quantified by 
the non-Gaussian parameter 



a 2 (t) 



3 ([r(f +t)-r(f )] 4 ) 
5 <[ r (f 0+ t)-r(i )] 2 } 2 



- 1. 



(12) 



Figure [5] shows a^if) for the oxygen atoms of PPO at var- 
ious temperatures. We see that the non-Gaussian param- 
eter exhibits a maximum in the late-/?/ early-a relaxation 
regime, consistent with findings for various model-glass 
formerSi 22 ' 31 i 32 ' 33 i 35 When the temperature is decreased, 
the position of the maximum, t a i, shifts to longer times. 
In Fig. [2ja) , we see that the temperature dependence of 
t a 2 is somewhat weaker, but still comparable to that of 
the a-relaxation time. The maximum 0:2(^02) increases 
upon cooling, in particular near the critical temperature 
T c , but the values are relatively small. More precisely, at 
comparable temperatures near T c , 0:2(^2) = 1-3— 2.5 was 
reported for the other studied models of glass-forming 
liquids, except for the silicon atoms in silica showing 
0:2(^2) = 0.8, see Ref. l35l for a detailed comparison. For 
PEO, the non-Gaussian parameter is even smaller than 
for PPO. This results in some ambiguities when extract- 
ing the maximum positions so that we refrain from dis- 
cussing the temperature dependence of t a i for PEO. 

To ascertain the spatial heterogeneity of the PEO and 
PPO dynamics, we demonstrate that highly mobile oxy- 
gen atoms form clusters larger than expected from ran- 
dom statistics. Following previous studieS ) 25 i 28 ' 29 i 33 ' 35 we 
characterize the particle mobility in a time interval t by 
the scalar displacement and select the 5% most mobile 
oxygen atoms for further analysis. Then, we define a 
cluster as a group of the most mobile oxygen atoms that 
reside in the first neighbor shells of each other. For both 
polymers, we use the position of the first minimum of 



FIG. 6: Mean cluster size S m (t) of highly mobile oxygen atoms 
for PEO at temperatures 450 K, 400 K, 350 K, 325 K, 300 K, 
and 280 K. When the particles used for the analysis are chosen 
irrespective of their mobilities, a mean cluster size of Sm « 4 
results, as indicated by the dashed line. 



the respective intermolecular oxygen-oxygen pair distri- 
bution function as criterion for the extension of the neigh- 
bor shell. Based on the probability distribution p s (n,t) 
of finding a cluster of size n for a time interval t, we 
calculate the weight-averaged mean cluster size 



S w (t) 



En n Ps( n .<) 



(13) 



This quantity measures the average size of a cluster to 
which one of the most mobile oxygen atoms belongs. Pre- 
viously, it was shown that the conclusions resulting from 
such analysis are unchanged when the fraction of highly 
mobile particles is varied in a meaningful range.— 

The mean cluster size S w (t) is displayed for PEO at 
various temperatures in Fig. [5] It is evident that S w (t) 
shows a maximum, which increases upon cooling. By 
contrast, a time and temperature independent small size 
S w f» 4 characterizes the clusters for the case of random 
statistics, i.e., when 5% of the oxygen atoms are chosen 
for analysis irrespective of their mobilities. These results 
clearly demonstrate the existence of spatially heteroge- 
neous dynamics. The transient nature of the clusters can 
be quantified, when we determine the times ts at which 
S w (t) is a maximum. In Fig. [51 we see for PEO and PPO 
that the clusters are largest in the a-relaxation regime 
at all studied temperatures. For PPO, we observe that 
the single maximum of S w (t) at high temperatures splits 
into two peaks near T C) see Fig. [7] This can be taken as 
another hint that a primary and a secondary relaxation 
coexist at sufficiently low temperatures. 

In Fig. [6l a closer inspection of the data for the highest 
temperatures reveals that S w (t) increases at long times 
in the diffusive regime. This increase is the mere con- 
sequence of the chain connectivity. At sufficiently long 
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FIG. 7: Mean string length L w (t) of highly mobile oxygen 
atoms for PPO at temperatures 450 K, 400 K, 350 K, 325 K, 
300 K, 280 K, 265 K, and 250 K. When the oxygen atoms used 
for the analysis are chosen irrespective of their mobilities, a 
mean string length of L w « 1.1 results, as indicated by the 
bottommost dashed line. The upmost dashed line is the mean 
cluster size Sw (t) of highly mobile oxygen atoms for PPO at 
T = 250K. 

times, the displacements of the individual atoms are 
largely determined by the displacement of the center-of- 
mass of the respective polymer chain and, hence, those 
oxygen atoms are highly mobile that belong to chains, 
showing the largest center-of-mass displacements from 
the statistical distribution. Due to their spatial proxim- 
ity along the chain, these atoms form extended clusters 
so that S w (t) increases when the single particle displace- 
ments start to become dominated by the center-of-mass 
displacements. By contrast, the maximum at shorter 
times is not due to chain connectivity. 

It has been shown that cooperative string-like mo- 
tion is an important channel for the relaxation of highly 
mobile particles in simple model-glass formersi 25 i 29 ' 30 i 31 
Therefore, we investigate the relevance of string-like mo- 
tion for the chemically realistic PEO and PPO models. 
Following these previous studies, we construct strings by 
connecting any two oxygen atoms i and j if the condition 

mm[\Ti(to)-Tj(to+t)\,\Tj(io)-Ti(to+t)\] < S 

holds for the atomic positions at two different times and 
set 5 to about 55% of the intermolecular oxygen-oxygen 
distance. Then, this condition means that one oxygen 
atom has moved and an another oxygen atom has oc- 
cupied its position. We checked that our conclusions are 
not altered when 5 is varied in a meaningful range. Using 
the above criterion, we determine the probability pi(l,t) 
of finding a string of length I for a time interval t and 
calculate the weight-averaged mean string length L w (t) 
in analogy with Eq. (fl~3]) . 

Figure[7]shows L w (t) for PPO at various temperatures. 
Evidently, the strings grow and shrink in time and they 
are substantially longer than that resulting from random 




FIG . 8 : Results for PEO atT = 280 K. (a) Typical trajectory 
<t>cc(t) of an OCCO dihedral angle during a time interval of 
1.15 ns. The straight lines mark the discontinuous trajectory 
s cc (t) resulting from mapping of the trajectory <f>cc(t) onto 
the conformational states g~ , t, and g + . (b) Probability dis- 
tributions of the OCCO (solid line) and COCC (dashed line) 
dihedral angles, p(<f> C c) and p(<j> oc ), respectively. 

statistics. Upon cooling, the position of the maximum, 
i^, shifts to longer times and the height of the maximum 
increases, where the growth is particularly prominent in 
the vicinity of T c . Thus, string- like motion is an impor- 
tant phenomenon at sufficiently low temperatures. L w (t) 
and S w (t) show a similar behavior. In particular, both 
quantities exhibit a two-peak signature at T = 250 K. The 
temperature dependence of tr, and t$ is compared in Fig. 
[21 For PEO and PPO, the strings, like the clusters, are 
largest in the a-relaxation regime. By contrast, simu- 
lation studies on simple glass-forming liquids, including 
a bead-spring polymer model, found that tr, and ts are 
located in the late-/3/ early-a relaxation regime, i.e., at 
significantly shorter times We conclude 

that the PEO and PPO models do show spatially het- 
erogenous and cooperative dynamics, however the char- 
acteristics of these effects differ from that for simple 
model-glass formers. 



D. Conformational relaxation 

Local conformational dynamics are of central impor- 
tance for the structural relaxation of polymers^ In the 
following, we investigate the torsional motion of the PEO 
model by analyzing the angular trajectories 4> C c(t) and 
</> oc (i), describing the time evolution of the OCCO and 
COCC dihedral angles, respectively. Figure [5] presents a 
typical trajectory 4> cc (t ). Comparison with the probabil- 
ity distribution of the dihedral angle, p (0 CC ), shows that 
the torsional motion is comprised of well defined transi- 
tions between the gauche^ trans (t), and gauche + 
(g + ) states, suggesting that a discretization is useful for 
an analysis of the conformational dynamics. Therefore, 
we map the trajectories 4> x (t) (x — cc,oc) onto discrete 
sequences s x (t ) of the conformational states, see Fig. [5] 
In this way, we eliminate effects from librational motions, 
which do not lead to structural relaxation. Then, analysis 
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FIG. 9: Time constants characterizing the conformational dy- 
namics of the OCCO dihedrals in the PEO model. We show 
the torsional correlation times t cc , defined as P cc (t cc ) — 1/e, 
and the rotational correlation times tr, see Fig. together 
with VFT interpolations. Moreover, the mean waiting times 
of the OCCO dihedrals in the t state are displayed. While all 
escape processes from the t state are considered for the calcu- 
lation of (t^ ), (fij?) ((tw b ) ) characterizes explicitly the forward 
(backward) jumps g^ 1 —yt—* (g ± ^t^g ± ). The solid lines 
are Arrhenius fits ({t* w ): £ a = 0.025 eV, {t*J): £ a = 0.044 eV, 
E a = 0.023 eV). Finally, we present the mean times (£3) 
needed for an OCCO dihedral to visit each of the three con- 
formational states at Jeast once. Calculated values of (£3) are 
shown as dotted and dashed lines. While g <-> g^ transitions 
were neglected when calculating the former data, these tran- 
sitions were taken into account when computing the latter 
data, see text for details. 



of the discrete sequences s x (t ) enables a straightforward 
characterization of the relevant torsional motion, e.g., in 
terms of waiting times t w and back-jump probabilities 
Pb, see below. 

To study the time scale of the conformational relax- 
ation, we use the discrete sequences s x and determine 
the probabilities p x (t) of finding a dihedral in the same 
conformational state at two times separated by a time 
interval t. Due to the finite number of conformational 
states, these probabilities exhibit a finite and temper- 
ature dependent plateau value p x (po). To remove this 
effect, we calculate the torsional correlation functions 



(Pxjt)} - (Pxjoo)) 
1 - (Px(00)) 



(14) 



and determine torsional correlation times t x according to 
Px(t x ) = 1/e. In Fig. [9j we compare the torsional corre- 
lation time T cc with the rotational correlation time tr, 
characterizing the reorientation of the C-H bond vectors. 
We see that both time constants show a comparable tem- 
perature dependence, implying that the conformational 
and the structural relaxations are related, consistent with 
results for other polymer models. 5,65 Some deviations in 
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FIG. 10: Probability distributions Pcc(tw) characterizing the 
waiting times of the OCCO dihedrals in the t state. Results 
for PEO at various temperatures are compared. 



the temperature dependence are expected since, in par- 
ticular at high temperatures, librational motions affect 
the rotational correlation times more than the torsional 
correlation times obtained from analysis of the discrete 
jump sequences. Moreover, the conformational dynam- 
ics of several dihedral species render the orientation of a 
given C-H bond time dependent and, hence, the absolute 
values of the rotational and torsional correlation times 
should differ. 

Next, we study the waiting times t w in the conforma- 
tional states, i.e., the time intervals between two subse- 
quent conformational transitions. We find that the wait- 
ing times differ between the OCCO and COCC dihedral 
angles and, for each dihedral species, they depend on 
the conformational state. Therefore, we separately de- 
termine the waiting times for each dihedral species and 
conformational state. The temperature dependent mean 
waiting time of the OCCO dihedrals in the t state, (i^), 
is included in Fig. [HI We see that the mean waiting time 
follows an Arrhenius law with a small activation energy 
E a = 0.025 eV. Qualitatively similar results are observed 
for all other dihedral species and conformational states. 
Hence, the temperature dependence of the mean waiting 
times is much weaker than that of the rotational and tor- 
sional correlation times, indicating that the longer resi- 
dence times in the conformational states at lower temper- 
atures are not sufficient to explain the strong slowdown 
of the a relaxation upon cooling, in harmony with results 
for other all-atom polymer models^ 

Closer insights into the nature of the conformational 
dynamics are available from the probability distributions 
of the waiting times, p x (t w ). Figure fTUl shows the distri- 
butions Pcciy w )i which characterize the waiting times of 
the OCCO dihedrals in the t state. In particular at the 
lower temperatures, there are substantial deviations from 
an exponential waiting-time distribution, indicating that 
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FIG. 11: Probability distributions p cc (i^) andpoc(im) charac- 
terizing the waiting times of the OCCO and COCC dihedrals 
in the t state, respectively. Results for PEO at T — 300 K are 
shown. For both dihedral species, we distinguish the waiting 



times prior to a forward jump (g - 
to a backward jump (g ± -^t~*g ± ). 



•t- 



j^) from that prior 



a Markov process does not apply to the conformational 
dynamics. At short waiting times, we see a nonexpo- 
nential decay, which exhibits a weak temperature depen- 
dence. At long waiting times, there is a nearly exponen- 
tial and more temperature dependent decay. Qualita- 
tively similar waiting-time distributions are observed for 
all dihedral species and conformational states, suggesting 
that the distributions are comprised of two contributions 
governing the behavior at short and long waiting times, 
respectively. 

In general, one can distinguish between forward and 
backward jumps. Let cri-i — ><7, — KTi+i be a sequence of 
three subsequently visited conformational states. Then, 
a forward and a backward jump in the state <jj are as- 
sociated with ^ai+i and o~i-\ =<7j_|_i, respectively. 
Unlike backward jumps, forward jumps enable visiting 
all three conformational states and, hence, they may be 
the cornerstone of conformational relaxation. Since tran- 
sitions <? ± ^.9 T are rare for the OCCO and COCC dihe- 
drals at the studied temperatures, the behavior in the t 
state is of particular importance for the exploration of all 
conformational states. Figure [IT] shows the probability 
distributions Px{t w ) characterizing the waiting times in 
the t state during the forward (g — >t^g T ) and back- 
ward (g — J >t^g ± ) sequences, respectively. We see that 
the results for the two dihedral species are comparable. 
Thus, it is not important that the t state is the majority 
state of 4>oc, whereas it is the minority state of </) cc , see 
Fig. [5] However, the waiting times strongly depend on 
the direction of the jumps. For both dihedral species, 
the forward jumps are characterized by nearly exponen- 
tial waiting-time distributions and, hence, the propensity 
to perform such transition does not depend on the jump 
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FIG. 12: Probability distribution p cc (iL) characterizing the 
waiting times of the OCCO dihedrals in the t state prior to a 
forward jump. Results for PEO at various temperatures are 
shown. The inset displays the temperature dependence of the 
probability pi that a forward jump g — > t is followed by a 
direct backward jump t—*g ± to the initial dihedral state. 



history. The backward jumps exhibit a more complex be- 
havior. In a semilogarithmic representation, the slope of 
the curves decreases with increasing waiting time until it 
becomes constant, indicative of an exponential behavior 
at sufficiently long waiting times. In the latter regime, 
the slope is comparable for the forward and backward 
jumps. 

These findings imply that both uncorrelated and cor- 
related conformational transitions occur. In some cases, 
the time and the direction of a jump are independent of 
the history, resulting in an exponential behavior and in 
similar rates of forward and backward jumps. In other 
cases, the dihedrals have a high tendency to return to the 
previous conformational state almost immediately, i.e., 
to perform a correlated backward jump, leading to the 
observed nonexponentiality at short waiting times. One 
can imagine that this proncness to a correlated back- 
ward jump depends on the jump history. For example, it 
may become weaker when time elapses after the previous 
transition. Also, a correlated backward jump may be less 
likely when a large number of forward-backward jumps 
have already taken place, i.e., when previous "unsuccess- 
ful" attempts have paved the way for a " successful" tran- 
sition to a new state, see below. 

To further investigate the role of the forward jumps 
for the conformational and structural relaxations, we an- 
alyze their temperature dependent behavior. Figure [T2] 
shows the distributions Pccitln) characterizing the wait- 
ing times of the OCCO dihedrals in the t state prior to 
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a forward jump at various temperatures. Evidently, all 
waiting-time distributions are nearly exponential, con- 
firming that the point in time of a forward jump is inde- 
pendent of the history. In Fig. we compare the tem- 
perature dependent mean waiting times (*£,'*) and (t^) 
characterizing the forward and the backward jumps in 
the t state, respectively. We see that the forward jumps 
exhibit a higher temperature dependence than the back- 
ward jumps. Specifically, the activation energies amount 
to E a = 0.044 eV for the former and E a = 0.023 eV for the 
latter. However, the temperature dependence of (tlf) is 
still much weaker than that of the torsional and rota- 
tional correlation times, corroborating our previous con- 
clusion that the longer waiting times at lower tempera- 
tures do not provide an explanation for the slowdown of 
the a relaxation upon cooling. 

These findings show that it is not sufficient to study 
single events, but it is necessary to correlate the times 
and the directions of two or more conformational transi- 
tions. First, we analyze the directions of two consecutive 
jumps. Specifically, we determine the back-jump proba- 
bility pi in the t state, i.e., the probability that a transi- 
tion 5 ± — »i is followed by a transition t—tg^. In Fig. [T2l 
the temperature dependent back-jump probabilities pi of 
the OCCO dihedrals are presented. We see a substantial 
increase of p b * upon cooling. Thus, when the tempera- 
ture is decreased, both longer waiting times and higher 
back-jump probabilities contribute to the slower explo- 
ration of the conformational states. In other words, the 
increase of the back-jump probabilities leads to an addi- 
tional delay of the conformational relaxation and, hence, 
to deviations from an Arrhcnius behavior. 

To study the exploration of the dihedral states in more 
detail, we determine the times tz needed for the dihedrals 
to visit all three conformational states. Specifically, we 
start from a randomly chosen time origin and define as 
the time elapsing until a given dihedral has first visited 
each of the three conformational states at least once. In 
Fig. [51 we include the mean times (t^) resulting from the 
conformational dynamics of the OCCO dihedrals at var- 
ious temperatures. Unlike the mean waiting times, the 
exploration times (£3) do not follow an Arrhcnius law. 
However, their temperature dependence is still weaker 
than that of the torsional correlation times, implying that 
the relation between the exploration process and the con- 
formational relaxation is not straightforward, but subtle. 
We note that previous work on polybutadiene did not 
relate the exploration process to the primary, but rather 
to the secondary relaxation; 7 

The question arises whether knowledge of the back- 
jump probabilities and of the mean waiting times prior 
to forward and backward jumps is sufficient to calculate 
the exploration times. To tackle this question, we assume 
that, at the randomly chosen time origin, an OCCO dihe- 
dral occupies one of the gauche states, being the majority 
states, and calculate the average time needed for the di- 
hedral to visit each conformational state at least once. 
First, we neglect direct transitions between the gauche 



states due to their rareness. Then, a forward jump in the 
t state is necessary to visit all dihedral states. Prior to a 
forward jump, the OCCO dihedral can perform back- 
ward jumps in the t state. Thus, in general, sequences 
g ± — >nbj(f— > g ± — ►) t^g T lead to the exploration of all 
conformational states. Provided the back-jump probabil- 
ity in the t state is independent of the history, the prob- 
ability of finding a sequence with ribj backward jumps in 
the t state is given by (1— pl)(pl) nbj . Then, the average 
time (tgtg) to move from one of the gauche states to the 
other can be written as 

(tgtg) = (1-P 6 *) E (Pb) nbj (tfl + n bj t b f) 

n b j =0 

= */< + t4-?^- (15) 

1 Pb 

Here, tbf is the time needed for a forward-backward 
jump sequence and tfi is the sum of the times elapsing 
prior to the first and the last conformational transition. 
When we assume that the waiting times can depend on 
the jump direction, but are otherwise independent of the 
history, tbf and tfi are determined by the mean waiting 
times (i^/) and (t^, h ) and by the mean waiting time (t 9 ^ ) 
in the g^ states, which is found to be essentially indepen- 
dent of the jump direction. Specifically, tbf=t^ J h +t 9 u and 
tfi = t 9 u /2+t^/ . The factor 1/2 in the latter equation is 
a consequence of the fact that, on average, the randomly 
chosen time origin lies in the middle of the waiting time 
before the first jump. 

Utilizing the knowledge of the mean waiting times and 
of the back-jump probabilities, we calculate the times 
(tgtg) according to Eq. (fT!)]) . In Fig. [SJ we see that the 
temperature dependence of (tgtg) deviates from an Ar- 
rhenius law due to the higher back-jump probabilities at 
lower temperatures. While the calculated and the actual 
exploration times are similar at high temperatures, the 
former show a weaker temperature dependence. There- 
fore, we dropped the assumption that direct transitions 
between the gauche states can be neglected and recalcu- 
lated the mean times (t g t g )- In Fig. [5J it is evident that 
considering the transitions between the gauche states has 
hardly any effect. Therefore, we refrain from specify- 
ing the equations for this case. Also, we expect that 
the assumption to start in one of the gauche states is 
not crucial, in particular at low temperatures, where the 
occupation of the trans state is small. The deviations 
between the calculated and the actual exploration times 
rather show that the assumption of history independent 
waiting times and back-jump probabilities is not justified, 
further illustrating the complexity of the conformational 
dynamics. 

Finally, we analyze the torsional motion of the OCCO 
dihedrals on a longer time scale t w 10 (t w ), where (t w ) 
is the mean waiting time resulting from all conforma- 
tional transitions of this dihedral species. We calculate 
the probability distribution z cc (n,t) of finding n confor- 
mational transitions during a time interval t. Figure [T^] 



11 




FIG. 13: Probability z cc (n,t) of Ending n conformational 
transitions of the OCCO dihedrals during a time interval t. 
We used t = 125.5 ps and t = 39.5 ps for PEO at T — 280 K 
({t w } = 14.3 ps) and T = 450 K ((t w ) = 3.9 ps), respectively. 
Thus, the ratio between the time interval and the mean wait- 
ing time amounts to t/t w ~ 10 for both temperatures. The 
solid line is the Poisson distribution calculated using the val- 
ues of t and t w at T = 450K, i.e., t/t w = 10.1. 



shows the distributions z cc (n,t) obtained for t « 10 (t w ) 
at T = 280 K and T = 450 K. The data for the higher 
temperature resemble the Poisson distribution for the 
used ratio tj (t w ) « 10 and, hence, a Markov process ap- 
proximates the conformational dynamics. By contrast, 
substantial deviations from the Poisson distribution are 
obvious for the lower temperature. Specifically, z cc (n,t) 
is much broader, indicating that a large fraction of di- 
hedrals performs less or more transitions than expected 
for the case of uncorrelated jump events. Thus, at suffi- 
ciently low temperatures, pronounced dynamical hetero- 
geneities do not only govern the structural relaxation, see 
Sec. IIII CI but also the conformational relaxation. 

Furthermore, it is evident from Fig. [T^] that it is more 
probable to find an even than an uneven number of con- 
formational transitions in a given time interval. This 
effect, which is most pronounced for small n and low 
temperatures, is another consequence of the existence of 
correlated forward-backward jumps. It implies that, at 
a given time, each OCCO dihedral has a preferred con- 
formational state. After the exit of its preferred state, a 
dihedral tends to return to this state within a very short 
period of time performing a correlated backward jump. 
Then, on average, the dihedrals spend much longer times 
in the preferred state than in the unpreferred state of the 
forward-backward sequence, resulting in a higher prob- 
ability of finding an even number of jumps, consistent 
with the observation. At T — 280 K, the difference be- 
tween even and uneven n is small for n > 30 and, hence, 
the information about the preferred state vanishes after 
about 30 conformational transitions, i.e., after about 15 
forward-backward jump sequences. 



E. Conclusions 

We have performed MD simulations to investigate the 
translational and rotational motion of chemically real- 
istic polymer models. For the studied PEO and PPO 
models, the structural relaxation shows the typical prop- 
erties of polymer melts, i.e., its time dependence dif- 
fers from a single exponential function and its temper- 
ature dependence deviates from an Arrhenius law, in 
harmony with results in the literature42 Specifically, for 
both polymer models, a KWW function well interpolates 
the translational and rotational correlation functions in 
the a-relaxation regime, time temperature superposition 
is obeyed, and a VFT law describes the temperature de- 
pendent a-relaxation times. The stretching and the tem- 
perature dependence of the correlation functions are con- 
sistent with results from experimental work, confirming 
the quality of the used force fields. 

Recently, the applicability of MCT to the structural re- 
laxation of chemically realistic polymer models was con- 
troversially discussed£ & 12 ' 13 Here, we have performed a 
MCT analysis for the PPO model using the incoherent 
intermediate scattering functions of the oxygen atoms. 
On the one hand, the analysis shows that a MCT power 
law with a critical temperature T c = 242 ± 2 K well de- 
scribes the temperature dependent correlation times in 
the range (T — T c )/T c = 0.1 — 0.9. Also, for a momentum 
transfer q corresponding to the intermolecular oxygen- 
oxygen distance, further MCT predictions are fulfilled. 
On the other hand, the temperature dependence of the 
a-relaxation time deviates from a MCT power law in 
the immediate vicinity of T c and the factorization theo- 
rem, being a central MCT prediction for the /3-relaxation 
regime, is violated. We conclude that MCT captures sev- 
eral aspects of the structural relaxation at appropriate 
temperatures and length scales, however it does not pro- 
vide a complete description since energy barriers against 
the torsional motion affect the dynamical behavior. 

Furthermore, we have studied the heterogeneity and 
the cooperativity of the structural relaxation. For the 
PEO and PPO models, we have demonstrated that highly 
mobile oxygen atoms aggregate into transient clusters, 
indicating the spatially heterogeneous nature of the dy- 
namics. Furthermore, we have shown that cooperative 
string-like motion facilitates the translational displace- 
ments of the highly mobile oxygen atoms at intermedi- 
ate times between ballistic motion and diffusive motion. 
Both clusters and strings increase in size upon cooling 
and, hence, spatial heterogeneity and cooperativity are 
prominent aspects of the molecular dynamics in partic- 
ular at low temperatures, in harmony with results for 
various models of glass-forming liquids^ Concerning the 
transient nature of these phenomena, we have found that 
the mean cluster size S w and the mean string length L w 
are a maximum at comparable times ts ~ tz in the a- 
relaxation regime for PEO and PPO. In this respect, the 
present findings differ from previous results. Specifically, 
the clusters and the strings were reported to be largest 
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at significantly earlier times in the late-/?/ early-a relax- 
ation regime for models of atomic liquids ; 25 ' 29 water—, 
silica ) 34 ' 35 and a bead-spring polymer . 28 ' 31 We conclude 
that the existence of spatially heterogeneous and cooper- 
ative dynamics are common to a broad variety of glass- 
forming liquids. However, the characteristics of these 
phenomena differ among the materials and, hence, relat- 
ing the properties of the heterogeneity and the cooper- 
ativity of the molecular dynamics to the respective fea- 
tures of the structural relaxation may yield interesting 
insights into the glass transition phenomenon. 

For the PEO model, we have demonstrated that a 
straightforward study of the conformational relaxation is 
possible, when mapping the continuous trajectories of the 
dihedral angles onto discrete sequences of the dihedral 
states. Various results have indicated a complex nature 
of the conformational dynamics. In particular, for both 
dihedral species and for all dihedral states, the probabil- 
ity distributions of the waiting times t w , i.e., of the time 
intervals between two subsequent conformational tran- 
sitions, strongly deviate from an exponential function, 
indicating that a Markov process does not apply to the 
conformational dynamics. To analyze the origin of this 
behavior, we have discriminated between backward and 
forward jumps, i.e., we have distinguished whether or not 
the dihedrals are in the same conformational state after 
exactly two transitions. This analysis revealed that the 
nonexponcntial waiting-time distributions are a conse- 
quence of correlated forward-backward jumps, which are 
an important aspect of the conformational dynamics at 
sufficiently low temperatures. 

Consistent with previous results for all-atom polymer 
models* 5 - we have observed that the mean waiting times 
show a weaker temperature dependence than the time 



constants of the conformational and structural relax- 
ations. Hence, it is not possible to explain the glassy 
slowdown on the basis of single events, but it is nec- 
essary to correlate the times and the directions of sev- 
eral consecutive conformational transitions. Analyzing 
the directions of subsequent transitions, we have shown 
that the probability of backward jumps increases upon 
cooling and, hence, the exploration of the conformational 
states is slower at lower temperatures due to both longer 
waiting times and higher back-jump probabilities. To 
obtain insights into the times of consecutive transitions, 
we considered probability distributions z x (n,t) of find- 
ing n transitions in a time interval t. In this way, we 
have revealed that the conformational dynamics resem- 
bles a Poisson process at high temperatures. However, 
when the temperature is decreased, dynamical hetero- 
geneities become important for the conformational re- 
laxation and the dihedrals start having a preferred con- 
formational state at a given time. After an exit of the 
preferred state state, the dihedrals show a high tendency 
to return to this state via a correlated backward jump. 
At the studied temperatures, the dihedrals remember this 
preferred state for up to about 30 conformational transi- 
tions. In other words, the conformational states are not 
sampled according to their statistical weights in the early 
stages of the conformational relaxation. 
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